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Fluctuations in the initial QCD matter density distribution are found to enhance the production 
of thermal photons significantly in the range 2 < pr < 4 GeV/c compared to a smooth initial state 
averaged profile in ideal hydrodynamic calculation for 200 AGeV Au+Au collisions at the Relativistic 
Heavy Ion Collider (RHIC) and 2.76 ATeV Pb+Pb collisions at the Large Hadron Collider (LHC). 
The thermal emission of photons is strongly dependent on the initial temperature of the system 
where the presence of 'hotspots' in the initial state translates into enhanced production of photons 
compared to a smooth profile. The effect of fluctuations in the initial state is found to be stronger 
for peripheral collisions and for lower beam energies. The pr spectra are found to be quite sensitive 
to the value of the initial formation time of the plasma which is not known unambiguously and 
which may vary with collision centralities at a particular beam energy. Increase in the value of the 
formation time lowers the production of thermal photons compared to the results from a shorter 
formation time. However, the relative enhancement from fluctuating initial states (compared to a 
smooth initial state) is found to be stronger for the larger values of formation time. The pt spectra 
alone are found to be insufficient to quantify the fluctuations in the initial density distribution due 
to the uncertainties in the initial conditions. A suitably normalized ratio of central-to-peripheral 
yield as a function of collision centrality and pr can be a useful measure of the fluctuation size scale. 



I. INTRODUCTION 

Ideal hydrodynamics with smooth initial density dis- 
tribution has been used successfully in recent past to 
model the evolution of the hot and dense system pro- 
duced in collisions of heavy nuclei at relativistic energies. 
The experimental data [l[ for px spectra and elliptic flow 
of several hadrons at the Relativistic Heavy Ion Collider 
(RHIC) energy are well reproduced upto a relatively large 
value of transverse momentum (pt ~ 2 GeV/c) using hy- 
drodynamic models 041] with smooth initial conditions 
(IC). However, for any given single event there may be 
randomly distributed inhomogenities in the initial energy 
density. Lately it has been shown that event-by-event 
(e-by-e) fluctuations in the IC have sizable effect on the 
physical observables if they are computed by averaging 
over many final states rather than considering an aver- 
aged initial density profile. 

Recent studies have reported that e-by-e fluctuating 
IC reproduces the experimental charged particle elliptic 
flow even for the most central collisions at RHIC which 
was underestimated by all earlier hydrodynamic calcula- 
tions using smooth IC. The e-by-e hydrodynamics with 
fluctuations in the IC also give a better agreement of the 
experimental charged particle spectra at high pt by in- 
creasing the number of particles there 0, H| • In addition, 
the initial state fluctuations may also help to understand 
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the various structures observed in two-particle correla- 
tions 0. 

Direct photons are considered as one of the few ob- 
servables able to probe the initial state of the hot and 
dense matter produced in the collisions. Thermal emis- 
sion of photons from the evolving system shows a very 
strong temperature dependence. Photons having pt > 
1 GeV/c are expected to provide a glimpse of the early 
part of the expansion history. These photons are es- 
pecially suitable for probing the fluctuations in the IC 
due to their sensitivity to the initial temperature of the 
plasma [Ioj |. 

Thermal radiation is predicted to be the dominant 
source of direct photons in the range 1 < px < 3 
GeV/c [ll|,[l^ ■ However, photon results using smooth IC 
and latest rates [HI, [l4| fall well below the experimental 
data points (Fig. 4 of Ref. (HI) in this pr range. In a re- 
cent work |l5( we have shown a substantial enhancement 
in the production of thermal photons for pj- > 1 GeV/c 
with an e-by-e fluctuating IC relative to a smooth initial- 
state averaged profile in the ideal hydrodynamic calcu- 
lation for 200 AGeV Au+Au collisions at RHIC and for 
the 0-20% centrality bin. Consequently, our results with 
fluctuating density profile improve the agreement with 
the PHENIX experimental data for 2 < pr < 4 GeV/c 
in that centrality bin. The enhancement was found to 
scale with the inverse of the initial fluctuation size. 

Here we perform a more systematic study of the ef- 
fects of the initial state fluctuations on the production 
of thermal photons in heavy ion collisions and propose 
an experimentally measurable quantity R2 p that can be 
useful to quantify the fluctuation size realized in nature. 

The effects of the initial density fluctuations are found 
to be more pronounced for lower beam energies and for 
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smaller size systems. In addition, the initial formation 
time of the plasma plays a crucial role in all these calcu- 
lations. As the value of To is not known unambiguously, 
we see that it is difficult to distinguish between the re- 
sults from fluctuating IC at a particular To and the re- 
sults at a relatively smaller To using the smooth IC. This 
is relevant since one can argue that the formation time of 
the plasma for peripheral collisions may not remain the 
same as for central collisions since the initial energy den- 
sity and system size become smaller moving from central 
to peripheral collisions. 

The effects of the initial state fluctuations on the pro- 
duction of thermal photons are studied with different val- 
ues of To changing with collision centrality at the LHC 
energy. The centrality dependence of To is taken from the 
EKRT model calculations [HI, [I?} ■ Larger values of To for 
peripheral collisions increase the difference between the 
results from the smooth and the fluctuating IC slightly 
compared to the results with a smaller To. From the cen- 
trality dependent results at different formation times we 
conclude that we need to look for a quantity besides the 
photon spectra as these are insufficient to quantify the 
fluctuations in the IC. The ratio of central to periph- 
eral yield normalized by the number of binary collisions, 
KJ p , shows potential to probe the density fluctuations 
and their size in the IC. 

The paper is organized as follows. First we briefly de- 
scribe the e-by-e hydrodynamic framework used in this 
study and the production of thermal photons from it. 
In the next section photon results are shown from the 
smooth and fluctuating IC at RHIC for different colli- 
sion centralities. After that we show the thermal photon 
Pt spectra from the smooth and fluctuating IC at the 
LHC energy and compare those with the RHIC results. 
Photon results from the fluctuating IC considering differ- 
ent values of the initial formation time To of the plasma 
at RHIC are shown in the next section along with the re- 
sults from centrality dependent To values at LHC. Next 
we define the ratio RJ p for thermal photons and calculate 
it as a function of collision centrality and the fluctuation 
size for different values of transverse momentum. Finally 
we summarize our results. 



II. THERMAL PHOTONS FROM 
EVENT-BY-EVENT HYDRODYNAMICS 

A. Event-by-event hydrodynamics 

We utilize here the e-by-e hydrodynamical model de- 
veloped in Q, which has been applied successfully to 
study the hadronic spectra and elliptic flow with fluc- 
tuating IC, and to calculate the thermal photon emission 
at RHIC ha and LHC energies. The model is based on 
ideal hydrodynamics with the assumption of longitudinal 
boost invariance and the remaining (2+l)-dimcnsional 
problem is solved numerically using the SHASTA algo- 
rithm [UBS. The E q uation of State (EoS), which shows 
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92 - 35 
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TABLE I: The A pa rt range, average impact parameter (&), av- 
erage numbers of participants and average numbers of binary 
collisions for different centrality classes of 200 AGeV Au+Au 
collisions at RHIC according to the MC Glauber model with 
&nn = 42 mb. 
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309.3 
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20 - 40 


220 - 104 


8.52 
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440.9 
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103 - 39 
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TABLE II: The iV P art range, average impact parameter (b), 
average numbers of participants and average numbers of bi- 
nary collisions for different centrality classes of 2.76 ATeV 
Pb+Pb collisions at LHC according to the MC Glauber model 
with cr at at = 64 mb. 



a sharp crossover transition from the plasma phase to 
hadronic matter, is taken from [20| . 

The IC are obtained from the Monte Carlo Glauber 
(MCG) model. The standard two-parameter Woods- 
Saxon nuclear density profile is used to randomly dis- 
tribute the nucleons in the colliding nuclei. Thus no 
nucleon-nucleon correlations or finite nucleon size effects 
are taken into account. Nucleons i and j from two differ- 
ent nuclei are assumed to collide whenever they satisfy 
the relation, 



(xi - Xjf + ( yi - Vj f < 



&NN 



7T 



(1) 



where <jnn is the inelastic nucleon-nucleon cross section 
and Xi , yi denote the positions of the nucleons in the 
transverse plane. The centrality classes arc defined by 
binning the distribution of events in the number of par- 
ticipants (A part ) as in Ref. Q- Table Q] shows the N pax t 
ranges corresponding to the different centrality classes, 
the average value of impact parameter and average num- 
ber of participants for each centrality class for 200 AGeV 
Au+Au collisions at RHIC. Table |H] shows the same 
for 2.76 ATeV Pb+Pb collisions at LHC. The inelastic 
nucleon-nucleon cross section is taken as 42 (64) mb for 
= 0.2 (2.76) TcV. 
The initial density profile is assumed to be proportional 
to the number of wounded nucleons (WN), where the 
entropy density s is distributed in the (jc, y) plane around 
the wounded nucleons using a 2D Gaussian, 
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We refer to this as the sWN profile. The parameter 
K is a fixed overall normalization constant and a is a 
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free parameter determining the size of the fluctuations. 
The effective interaction radius of the colliding nucleons, 
\Z&nn /tt /2 ~ 0.6 fm (for RHIC), sets a natural order of 
magnitude for the size parameter. The default value of 
a is taken as 0.4 fm from Ref. @, [Hj], however, we vary 
the value of a from 0.4 to 0.8 fm to check the sensitivity 
of the results to the fluctuation size. 

Both intra-event fluctuations (i.e., 'hotspots' and 
'holes') and the inter-event fluctuations in the total en- 
tropy between different events arc thus generated by the 
MCG modeling of the initial state. We call events with 
larger than average entropy "hot events" and those with 
smaller than average entropy " cold events" . 

The constant K is taken as 102.1 fm -1 for Au+Au 
collisions at ^/snn = 200 GcV. First the initial formation 
time (tq) of the plasma is fixed at To = 0.17 fm/c (as in 
Refs. M.IT5I0 . motivated by the EKRT minijet saturation 
model [16| for all the events. However, in the later part of 
this study we vary the value of tq with collision centrality 
to check the effect of these variations and fluctuations 
in the IC for thermal photon radiation. The freeze-out 
temperature is 160 Me V, which reproduces the measured 
Pt spectrum of pions [2l[ well for the sWN profile for a 
fixed To. 

For the 2.76 ATeV Pb+Pb collisions at LHC we use the 
initial formation time To = 0.14 fm/c taken from EKRT 
model (l7| . The overall normalization constant in this 
case is K = 250 fm -1 which reproduces the measured 
multiplicity of about 1600 (22[. The freeze-out tempera- 
ture is kept same at RHIC and LHC. 



B. Thermal photon emission 

The quark-gluon Compton scattering and quark-anti- 
quark annihilation are the leading order processes for 
thermal photon production in the partonic phase. A 
significant contribution also comes from the resummed 
next to leading order bremsstrahlung processes [HI . The 
rate of thermal photon production in the QGP depends 
on the momentum distribution of the partons which are 
governed by the thermodynamical conditions of the mat- 
ter. 

The 7r and p mesons contribute dominantly to the pho- 
ton production from a hot hadronic gas having temper- 
ature of the order of pion mass. This is due to the low 
mass of pions and the large spin iso-spin degeneracy of p 
mesons, makin g th em the most easily accessible particles 
in the medium [23J . The leading photon producing chan- 
nels involving 7r and p mesons are tttt — > pr/, irp — > nj, 
and p — > 77777. 

We use the rates R = EdN/d 3 pd A x of for the 
plasma and those of [l4| for the hadronic matter which 
at present can be considered as the state of the art. The 
transition from the plasma rates to the hadronic rates is 
assumed to happen at a temperature of 170 MeV in this 
study. The total thermal emission from the quark and 
hadronic matter phases is obtained integrating the rates 



over the space-time evolution of the medium, 

EdN/d 3 p = J d i xR^E*(x),T{x)y (3) 

where E*{x) = p fl u fl (x). The 4-momentum of the pho- 
ton is p* 1 = {pt cosh Y,pt cos</>, pr sin <fi,pT smhY), and 
the 4- velocity of the flow field is u M = 7t (cosh 77, v x , v y , 
sinhri) with jt = (1— ^t)" 1 ^ 2 , v t — v x+ v y- The vol- 
ume element is d A x = t dr dx dy drj, where t = (i 2 — z 2 ) 1 / 2 
is the longitudinal proper time and 77= tanh -1 (z/t) is 
the space-time rapidity. The photon momentum is 
parametrized by its rapidity Y", transverse momentum 
Pt, and azimuthal emission angle <p. 

C. Averaging procedure 

The smooth initial density distribution from the MC 
Glauber calculation is obtained by taking an average of 
1000 fluctuating initial states The sWN profile 

must be converted to energy density distribution using 
the EoS, since the input to hydrodynamics is the energy 
density profile. When averaging over events in the ini- 
tial state, we apply the EoS in each event and in the 
end the average is taken over energy densities. Although 
photons are emitted during the entire space-time evolu- 
tion, the yield for pt > 1 GeV/c is dominated by QGP 
radiation. Thus averaging over 1000 events is enough to 
remove the fluctuations from the interior region of the 
fireball. 

The results from the fluctuating IC are obtained by 
averaging the photon spectra over a sufficiently large 
number of random events where the number needed in- 
creases for peripheral collisions. For example, we choose 
40 events for 0-20% centrality bin and 80 events for 40- 
60% centrality bin. In order to fix the necessary number 
for all the centrality bins we check that addition of an- 
other very hot or very cold event does not change the 
results significantly. 

III. RESULTS FROM FLUCTUATING AND 
SMOOTH IC 

A. Centrality dependence at RHIC 

First we study the centrality dependent results from 
200 AGcV Au+Au collisions at RHIC considering To = 
0.17 fm/c. Figure Q] shows the comparison between the 
smooth and the fluctuating IC for thermal photon px 
spectra at different collision centralities with a = 0.4 fm. 
The solid lines with closed symbols in the figure show 
the results from the fluctuating IC and the dashed lines 
with open symbols are the results from a smooth initial 
state averaged IC (See Table [T] for the A" part ranges and 
the corresponding values of average impact parameter in 
the different centrality bins). 
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FIG. 1: (Color online) Centrality dependence of thermal pho- 
ton pr spectra for 200 AGeV Au+Au collisions at RHIC from 
fluctuating and smooth IC. 




p T (GeV/c) 



FIG. 2: (Color online) Thermal photon pr spectra from 
smooth and fluctuating IC for 2.76 ATeV Pb+Pb collisions 
at LHC at different collision centralities and at size parame- 
ter a — 0.4 fm. 



We see from Fig. [T] that the enhancement in the pro- 
duction from the fluctuating IC grows for peripheral col- 
lisions compared to the results from the smooth IC. The 
exponential slope of the spectra from the fluctuating IC 
is found to be about 10%, 12% and 16% flatter than what 
is obtained from the smooth IC for - 20 %, 20 - 40 % 
and 40 - 60% centrality bins respectively, in the region 
2 < pt < 4 GeV/c. Moving from central to peripheral 
collisions, the total number of 'hotspots' decreases in the 
initial density distribution as the number of participants 
is less for peripheral collisions. However, the presence 
of even a few 'hotspots' in initial state makes significant 
enhancement in the production for peripheral collisions 
as the photon yield is very small there compared to the 
central collisions. Thus the relative importance of initial 
state fluctuations increases towards peripheral collisions. 

We have shown in the earlier work [l5j that production 
of thermal photons shows a very strong dependence on 
the initial state density fluctuations. The 'hotspots' in 
the fluctuating IC enhance the production significantly 
in the range pt > 1 GeV/c compared to the smooth ini- 
tial state averaged IC for 0-20% centrality bin at RHIC. 
This is due to the exponential temperature dependence of 
photon emission; the high pt photons arc mostly emitted 
from the very early stage of the system expansion when 
hydrodynamic flow is still weak. The low px photons 
are mainly from the relatively cold later stage of the sys- 
tem expansion and are not effected significantly by the 
presence of 'hotspots' in the IC, as the strong pressure 
gradients remove the 'hotspots' rapidly. 

Wc have checked that lowering the value of the freeze- 
out temperature from 160 to 120 MeV enhances the pro- 
duction significantly in the hadronic phase. This addi- 
tional contribution to the total spectrum, however, is 
not substantial for pr > 2 GeV/c for central collisions. 
The relative contribution from the hadronic phase com- 
pared to the QGP phase increases towards peripheral col- 
lisions. Thus, the lowering of the freeze-out temperature 
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FIG. 3: (Color online) Ratio of the thermal photon spectra 
from the smooth and the fluctuating IC for 0-20% centrality 
bins at RHIC and LHC energies and for size parameter a = 
0.4 fm. 



increases the production from hadronic phase in the low 
Pt region more for peripheral collisions. However, fluc- 
tuations in the IC enhance the production even more 
for peripheral collisions and a smaller value of freeze-out 
temperature does not change results significantly in the 
thermal px range. It is to be noted that a bag model 
EoS produces more hadronic photons near the transition 
region compared to a lattice based EOS, however the 
contribution from the hadronic phase is not si gnif icant 



for p T > 2 GeV/c even for a bag model EOS [24j. In 
addition, wc found that these results are quite sensitive 
to the initial profile details. The sWN profile leads to 
more pronounced maxima and hotter initial state than 
the energy initialization (eWN) profile [l5|. Thus also 
the thermal photon pt spectra are harder with the sWN 
profile. 
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FIG. 4: (Color online) Thermal photon pr spectra at differ- 
ent centrality bins for to values of 0.17 fm/c (gray lines) and 
0.2 fm/c and at a fixed value of size parameter a = 0.4 fm. 
Results from fluctuating and smooth IC are shown by solid 
lines with closed symbols and dashed lines with open symbols 
respectively. 



B. Results at LHC and comparison with RHIC 
results 

From the RHIC results at different collision centralitics 
we saw that the effect of the initial state fluctuations is 
stronger when the system is smaller in size and the initial 
energy is less. Thus, one can expect that for Pb+Pb 
collisions at the LHC, where the initial energy density 
as well as the temperature is significantly higher than at 
RHIC, the effect of fluctuations will be less pronounced. 

Figure [2] shows the thermal photon spectra for different 
collision centralities at the LHC (See Table |H] for the 
^Vpart ranges). Results from the smooth IC arc shown 
by the dashed lines with open symbols and those from 
the fluctuating IC by the solid lines with closed symbols. 
The production of thermal photons is found to increase 
by about a factor of 4 at pt =1 GeV/c to a factor of 
more than 10 at p T =4 GeV/c from RHIC to LHC for 
the 0-20% centrality bin. The difference between the 
fluctuating and smooth initial conditions shows a similar 
centrality dependence at LHC and at RHIC. 

The ratio of thermal photon yield from fluctuating and 
smooth IC for 0-20% centrality bin at the LHC is shown 
in the Fig. [3] where the result is compared with the RHIC 
result from Au+Au collisions at the same centrality bin. 
For both the energies the value of the size parameter a 
is taken as 0.4 fm. As seen in Figs. [T] and f2J fluctua- 
tions in the initial state enhance the production more for 
peripheral collisions than for central collisions. However, 
the overall magnitude of the enhancement is less at the 
LHC than at RHIC in the same pr range. 

The results at LHC along with the centrality depen- 
dence results at RHIC confirm the earlier expectation 
that one might have to look at lower beam energies and 
smaller size systems for more sensitive measurements of 
the initial state fluctuations in heavy ion collisions. 



C. Formation time dependence of the results at 
RHIC 

The initial formation time plays a very crucial role in 
the production of thermal photons from the plasma phase 
where a smaller tq leads to a larger initial temperature 
and enhanced production of high px photons [2 51 ] . Thus, 
a significant enhancement in the photon production can 
also be obtained by reducing the value of the initial for- 
mation time for the smooth IC. As the value of the initial 
formation time is not known precisely, it is important to 
know whether the enhancement of the photon spectra 
due to fluctuations in the IC can be distinguished from 
the enhancement due to a smaller value of the formation 
time. A systematic study of thermal photon spectra at 
different centrality bins and at different values of tq is 
therefore necessary. 

Figure |4] shows the thermal photon spectra from the 
fluctuating and smooth IC at RHIC for two different val- 
ues of To; 0.17 and 0.20 fm/c, where the value of r is 
changed by keeping the total entropy of the system fixed. 
We see that even a small change in the value of the for- 
mation time from 0.17 to 0.20 fm/c affects the results 
from both the smooth and the fluctuating IC visibly. In 
the fluctuating case, the effect is stronger for the central 
collisions than for peripheral collisions as the relative en- 
hancement in the production due to fluctuations is less 
for central collisions. 

The value of to differs significantly in different hydro- 
dynamic models 0, @| and mostly it varies from 0.17 to 
0.60 fm/c for Au+Au collisions at maximum RHIC en- 
ergy. We consider here 0.17 and 0.60 fm/c as the lower 
and upper limits of tq at RHIC and compare the photon 
spectra computed with these values. 

When the initial time To or the fluctuation size param- 
eter <j is changed, also the freeze-out temperature should 
be tuned in order to fit the hadron spectra. Here we how- 
ever keep the freeze-out temperature fixed for simplicity 
since the photon emission from later stages of hadron 
gas phase is so small that total yield in the px range con- 
sidered here would not change significantly. Note that 
this might not be true for other variables such as photon 
elliptic flow 

The size parameter a is a free parameter in our calcu- 
lation and it is thus interesting to the see the effect of 
a simultaneous change in the formation time and in the 
value of a on the production of thermal photons from 
the fluctuating IC. The upper panel of Fig. O shows the 
photon pt spectra from the fluctuating and smooth IC at 
To values of 0.17 (gray lines) and 0.6 fm/c (colored lines) 
for different collision centralities and at a fixed value of 
the size parameter a = 0.4 fm. The lower panel of the 
same figure shows the results for the same parameters 
but with <7 = 0.8 fm. The yield falls sharply in the range 
2 < pt < 4 GeV / c for all collision centralities as well as 
for both the a values when To is increased from 0.17 to 
0.60 fm/c. For a particular centrality bin, the difference 
between the smooth and the fluctuating IC is increased 
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FIG. 5: (Color online) Thermal photon pr spectra at different 
centrality bins for to values of 0.17 fm/c (gray lines) and 0.6 
fm/c (colored lines) and for size parameter a — 0.4 (upper 
panel) and 0.8 (lower panel) fm. Results from fluctuating 
and smooth IC are shown by solid lines with closed symbols 
and dashed lines with open symbols respectively. 



for larger values of to . This is again due to the fact that a 
larger value of t leads to a relatively smaller value of ini- 
tial temperature which is similar to a peripheral collision 
with smaller value of To. 

These results show that it is difficult to conclude any- 
thing about the fluctuations in the IC and the size pa- 
rameter a as long as tq is unknown. 



IV. RESULTS WITH CENTRALITY 
DEPENDENT t 

So far we have considered different collision ccntralitics 
by keeping the initial time t fixed. However, in periph- 
eral collisions the system is more dilute and the average 
energy scale of the produced quanta can be expected to 
be smaller than in central collisions, and thus both the 
system's formation and thermalization should take longer 
than in central collisions. Consequently, one should con- 
sider a To which is systematically increasing from central 
to peripheral collisions [26|, [27|. With the same logics, 
even within one collision, the tq should change according 



to the initial density fluctuations in the system. Clearly, 
such an advanced modeling of a fully local To is not yet 
available. 

To get an order of magnitude estimate for the changes 
in the average t from one centrality bin to another, we 
may, nevertheless perform an LHC case-study by using a 
simple non-local version of the EKRT model Ml: First, 
as in [H, [28[ , the average number of participants for each 
centrality bin is obtained from the optical Glauber model, 
and an effective mass number A e g < A is determined so 
that the number of participants in a central A c ff+A c g 
collision matches with that of the centrality class in ques- 
tion. Then, an average To for each centrality class is ob- 
tained by requiring that the production of minijets (glu- 
ons) of transverse momenta px > po 3> Aqcd in the 
mid-rapidity unit of the central A c fi+A jf collision fill 
the available effective transverse area, 



N Aclt A ott (b = 0,p , Ay = 1, Vs) ■ ir/pl 



(4) 



where Ra^ is the radius of the effective nucleus, Na b{1 a b{1 
is the number of minijets computed in LO pQCD 1 , and 
tt/pq is the effective uncertainty area assigned to each 
minijet. The above condition then gives the saturation 
momentum as po = p sa t and the average formation time 
of the system as l/p sa t- Assuming an instant thermal- 
ization - which we believe is a reasonable approximation 
since the build-up of pressure and flow in the system 
starts very early - one arrives at the To values given in 
Table III 2 . For clarity, we note that the multiplicities at 
different centralities are still fixed here as explained in 
Sees. II. A. and III. C; a more systematic study of mul- 
tiplicities in the NLO EKRT framework will be presented 
elsewhere [l7| . 

Tabic |lll] shows the obtained values of To, Asff and p sa t 
for three different centrality classes at LHC. We notice 
that the change in To from the 0-5 % to the 40-50 % cen- 
trality bin is only of the order of 30 %, i.e. much less than 
the overall uncertainty (factor ~ 3.5) in the To considered 
in the previous section. Figure [5] shows the corresponding 
thermal photon px spectra at the LHC from the smooth 
and fluctuating IC for 20-30 %, and 40-50 % centrality 
bins using the centrality-dependent t along with the re- 
sults obtained for a constant To = 0.14 fm/c (as obtained 
for central collisions here). As we can see in the figure, 
for peripheral collisions and fixed To the fluctuations in 
the IC enhance the production significantly (from the 
gray dashed lines to the gray solid lines). This enhance- 
ment is, however, partially compensated by a suppression 
caused by a larger formation time (from the gray solid 
lines to the colored solid lines). As a net result a visi- 
ble enhancement relative to the case with smooth IC and 



1 Without any ^-factors or an EKRT-based fit to the LHC mul- 
tiplicity |29|| . as we are only after a rough to estimate here. 

2 We thank Risto Paatelainen for computing these numbers for us 
from a work in progress [l7ll . 
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FIG. 6: (Color online) Thermal photon spectra for periph- 
eral collisions at the LHC from fluctuating and smooth IC 
considering centrality dependent to from the EKRT model. 



Centrality bin (%) 


Po (GeV) 


A cfi 


r (fm) 


0-5 


1.39 


193 


0.140 


20 - 30 


1.21 


90 


0.164 


40 - 50 


1.05 


40 


0.188 



TABLE III: value of to for there different centrality bins 
for 2.76 ATeV Pb+Pb collisions at LHC calculated from the 
EKRT model (see Ref. p]J for details). 



fixed To (from the dashed gray to the solid colored lines) 
remains. 

From all the results at RHIC and LHC above, one can 
conclude that it is difficult to distinguish between the ef- 
fects of the IC fluctuations and formation time only by 
looking at the pr spectra of thermal photons. Instead, 
one needs to look for some other experimentally measur- 
able quantity where the ambiguities regarding the initial 
conditions can be reduced. 



Centrality % 


Apart range 


(b) (fm) 


(Apart) 


(iVbin) 


- to 


394 - 276 


3.14 


326.3 


958.9 


10 - 20 


275 - 197 


5.72 


234.0 


599.0 


20 - 30 


196 - 138 


7.40 


165.6 


369.6 


30 - 40 


137 - 93 


8.76 


113.9 


217.8 


40 - 50 


92 - 59 


9.94 


74.55 


119.9 


50 - 60 


58 - 34 


11.0 


45.20 


59.74 



TABLE IV: -/Vpart, average impact parameter, average number 
of participants and average number of binary collisions for 
different centrality bins and for 200 AGeV Au+Au collisions 
at RHIC with a NN = 42 mb. 
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FIG. 7: (Color online) R2 P as a function of collision central- 
ity for thermal photons at different pr values for 200 AGeV 
Au+Au collisions at RHIC. Two different values of a, 0.4 and 
0.8 fm are used. The ./Vp ar t ranges for the different centrality 
bins are taken from Tabic ITvl 



A. Ratio of central to peripheral yield at RHIC 

The ratio of direct photon spectra from central to 
peripheral collisions is a useful measure since all overall 
normalization factors (which may not be well known) are 
cancelled. This ratio as a function of collision centrality 
also offers the possibility to better focus on the system 
size dependence of the IC fluctuation. 

We define a quantity RJ p , which is the ratio of central 
to peripheral yield of thermal photons normalized by the 
number of binary collisions 



R1 | = dN/d 2 p T dY\ _ 10% ^ N bin \i-j% 
cp 1 dN/d 2 pTdY\i_ j% N bi n\o-w% ' 



(5) 



where the value of i is changed from 10 to 70 in steps of 
10 and j=i+10 in the following. Note that although not 
considered here, for pQCD direct photons such ratio is 
close to unity. We calculate RJ p as a function of collision 
centrality for different values of transverse momentum as 
well as a function of the size parameter a. 

Note that we choose a different definition of RJ„ rather 

IT 

than the conventional definition where the result from 
the most peripheral collision is kept as a fixed denomina- 
tor and the numerator is changed for different centrality 
bins. The hydrodynamical framework is more realistic 
for central collisions than for peripheral collisions and 
as a result we keep the numerator fixed and change the 
denominator accordingly. 

Figure [7] shows the variation of RJ p as a function of 
collision centrality for different values of transverse mo- 
mentum and for a values of 0.4 and 0.8 fm. The results 
both from the smooth and the fluctuating IC approach 1 
towards central collision by definition and for a fixed px a 
clear difference between the results from the fluctuating 
and smooth IC can be observed for peripheral collisions. 

At pt = 1 GeV/c, differences between the results from 
the fluctuating and the smooth IC are small for central 
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Centrality % 


iVpart range 


(b) (fm) 


(Apart) 


(iVbin) 


0-10 


416 - 306 


3.27 


357.5 


1506 


10 - 20 


305 - 221 


6.00 


261.0 


929.0 


20 - 30 


220 - 155 


7.78 


186.1 


562.3 


30 - 40 


154 - 104 


9.26 


127.7 


319.1 


40 - 50 


103 - 66 


10.5 


83.36 


168.8 


50 - 60 


65 - 39 


11.6 


51.34 


82.37 



TABLE V: iVp ar t, average impact parameter, average number 
of participants and average number of binary collisions for 
different centrality bins and for 2.76 ATeV Pb+Pb collisions 
at LHC with onn = 64 mb. 



10 20 30 40 50 10 20 30 40 50 60 
centrality (%) 



FIG. 8: (Color online) R2 P as a function of collision central- 
ity for thermal photons at different px values for 200 AGeV 
Au+Au collisions at RHIC. Size parameter a is kept fixed at 
0.4 fm. 
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FIG. 9: (Color online) R2 P at px values of 2, 3, 4, and 5 GeV/c 
for 200 AGeV Au+Au collisions at RHIC as a function of size 
parameter a from fluctuating IC. 



collisions and they increase slightly for peripheral colli- 
sions. In the thermal region 2 < px < 4 GeV/c, RJ p from 
the fluctuating IC does not change significantly with colli- 
sion centrality. In contrast, for smooth IC where only the 
To dependence is relevant, the value of RJ p rises rapidly 
with collision centrality for both values of a. Fluctua- 
tions in the IC enhance the production more for periph- 
eral collisions compared to the smooth profile case. As a 
result the central to peripheral ratio for the smooth IC 
is always larger than for the fluctuating IC and the dif- 
ference increases towards the larger values of transverse 
momentum. 

For larger values of a, the effect of fluctuations becomes 
less pronounced and as a result RJ p for peripheral colli- 
sions increases with the centrality at large px- Results 
from the smooth IC rise faster compared to the fluctu- 
ating IC at a — 0.8 fm indicating that some fluctuation 
effects are still present even with this a value. 



In Fig. [5J RJ p for thermal photons as a function of colli- 
sion centrality is shown for To = 0.6 fm/c and at different 
values of px- For the smooth IC, the value of RJ p in- 
creases rapidly from central towards peripheral collisions 
and the increase is more for larger values of transverse 
momentum, whereas for the fluctuating IC the variation 
of R2 p with collision centrality is not very rapid, and only 
for a very large px a significant increase in the value of 
RJ p can be observed. Thus, comparing Figs. [7] and [5J it 
can be seen that the difference between the RJ„ values 
from the smooth and fluctuating profiles is more pro- 
nounced for larger tq. These results show that R1 p is a 
potentially useful observable in the study of IC fluctua- 
tions and the To uncertainty. 

The RJ„ results for four different values of transverse 
momentum are shown in Fig. [9] as a function of the size 
parameter a at different collision ccntralitics. For px =2 
GeV/c RJ p shows very little variation with a for almost 
all the centralities as the 'hotspots' in the fluctuating IC 
affect the spectra marginally in this px range. However 
for px = 3 GeV/c, RJ p shows stronger dependence on a 
as fluctuations in the IC affect the spectra more signifi- 
cantly and the effect is stronger for smaller values of a. 
A much larger variation of RJ p with a can be seen for 
px > 4 GeV/c and even for central and mid-central col- 
lisions. We have good reasons to believe that the direct 
photon spectra in the px range 3-4 GeV/c is dominated 
by thermal radiation [151 | and the photon spectra can be 
measured experimentally up to a very large collision cen- 
trality as shown in the direct photon data by PHENIX. 
Thus, an experimental measurement of the parameter 
RJ p at different px values and at different collision cen- 
tralities can actually serve as a probe for the size param- 
eter a. 

Finally, in Fig. [TU]we plot the central to peripheral ra- 
tio at LHC (see Table |V] for the 7V part and iVbi n ranges for 
different centrality bins at LHC) . We notice that the ratio 
R2 p is generally much flatter at higher collision energy. 
Also the difference between the smooth and fluctuating 
initial conditions is smaller at the LHC than at RHIC. 
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FIG. 10: (Color online) R2 P as a function of collision central- 
ity for thermal photons at different pr values for 2.76 ATeV 
Pb+Pb collisions at LHC. Two different values of a, 0.4 and 
0.8 fm are used. N pa ,it ranges are taken from Table IVl 

V. SUMMARY AND CONCLUSIONS 

In conclusion, we have shown that fluctuations in the 
IC lead to a significant enhancement in the production 
of thermal photons compared to a smooth initial state 
averaged IC in the ideal hydrodynamic calculation. For 
a particular beam energy the enhancement is found to be 
more pronounced for peripheral collisions than for cen- 
tral collisions. The relative enhancement is found to be 
comparatively less at LHC than at RHIC for the same 
centrality bin. 

The thermal photon pt results at RHIC and LHC are 
found to be quite sensitive to the value of the initial for- 
mation time. A larger value of formation time lowers the 
production of thermal photons generically. One might 
also expect that the value of To is not fixed and it may 
change with collision centralities. We have shown the pr 
spectra from the smooth and fluctuating IC considering 
different values of the initial formation time as well as re- 



sults considering ccntrality dependent tq values from the 
EKRT model. The effect of fluctuations in the IC and a 
larger value of formation time act in opposite directions 
by increasing and decreasing the production respectively. 
All these results suggest that the effect of the IC fluctu- 
ations and change in the formation time can not be dis- 
tinguished easily only by studying the spectra of thermal 
photons. 

We have shown that the central to peripheral ratio of 
Pt spectra (normalized with the number of binary col- 
lisions) can be a useful measure to reduce the model- 
dependent uncertainties in the calculations. Experimen- 
tal determination of RJ p at different pt and at different 
collision centralities would be a valuable tool for esti- 
mating the size scale of the fluctuations. The difference 
between the RJ n results from the smooth and the flue- 
tuating IC has been found stronger for peripheral colli- 
sions and for larger values of pr- We also found that at 
the LHC energy the RJ p is much flatter than at RHIC 
and the difference between fluctuating and smooth IC 
is smaller. However, one should bare in mind that hy- 
drodynamics is not expected to work at very peripheral 
collisions (> 60%), and that at large pr (> 4 GeV/c) the 
pQCD photons dominate the direct photon spectra. 



Acknowledgments 

We thank R. Paatelainen for providing us with his 
EKRT code and the To values for peripheral collisions at 
LHC. We gratefully acknowledge the financial support by 
the Academy of Finland. TR and RC are supported by 
the Academy researcher program (project 130472) and 
KJE by a research grant (project 133005). In addition, 
HH was supported by the national Graduate School of 
Particle and Nuclear Physics and the Extreme Matter 
Institute (EMMI). We acknowledge CSC - IT Center for 
Science in Espoo, Finland, for the allocation of compu- 
tational resources. 



[1] C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. [6] 
87, 182301 (2001); ibid 89, 132301 (2002); ibid 90, 
032301 (2003); S. S. Adler et al. [PHENIX Collabora- 
tion], Phys. Rev. Lett. 91, 182301 (2003). [7] 

[2] P. Huovinen, P. Kolb, U. Heinz, P. V. Ruuskanen, and 

S. Voloshin, Phys. Lett. B 503 , 58 (2001); D. Teaney, [8] 
J. Lauret, and E. Shuryak, |rmcl-th/0110037| 

[3] H. Niemi, K. J. Eskola, and P. V. Ruuskanen, Phys. Rev. [9] 
C79, 024903 (2009). 

[4] K. J. Eskola, H. Honkanen, H. Niemi, P. V. Ruuska- 
nen and S. S. Rasanen, Phys. Rev. C 72, 044904 (2005); 
P. Huovinen and P. V. Ruuskanen, Ann. Rev. Nucl. Part. [10] 
Sci. 56, 163 (2006); C. Nonakaand S. A. Bass, Phys. Rev. 
C 75, 014902 (2007). [11] 

[5] T. Hirano, Phys. Rev. C 65, 011901 (2002); T. Hirano [12] 
and K. Tsuda, Phys. Rev. C 66, 054905 (2002). 



P. F. Kolb and U. W. Heinz, Quark gluon plasma 3, ed. 
R. C. Hwa and X. N. Wang (Singapore: World Scientific) 
634 (2003) [nucl-th/0305084] . 

H. Holopainen, H. Niemi, and K. J. Eskola, Phys. Rev. 
C 83, 034901 (2011). 

Y. Hama, T. Kodama, and O. Socolowski, Jr., Braz. J. 
Phys. 35, 24 (2005). 

R. Andrade, F. Grassi, Y. Hama, T. Kodama, and 
O. Socolowski, Jr., Phys. Rev. Lett. 97, 202302 (2006); 
R. P. G. Andrade, F. Grassi, Y. Hama, T. Kodama, and 
W. L. Qian, Phys. Rev. Lett. 101, 112301 (2008). 
P. V. Ruuskanen, Nucl. Phys. A 544, 169 (1992), and 
references therein. 

D. K. Srivastava, J. Phys. G 35, 104026 (2008). 

A. Adare et al. [PHENIX Collaboration], Phys. Rev. 

Lett. 104, 132301 (2010). 



10 



[13] P. Arnold, G. D. Moore, and L. G. Yaffe, JHEP 0112, 
009 (2001). 

[14] S. Turbide, R. Rapp, and C. Gale, Phys. Rev. C 69, 
014903 (2004). 

[15] R. Chatterjee, H. Holopainen, T. Renk, and K. J. Es- 
kola, Phys. Rev. C83, 054908 (2011); R. Chatterjee, 
H. Holopainen, T. Renk, K. J. Eskola, J. Phys. G 38 
124136 (2011). 

[16] K. J. Eskola, K. Kajantie, P. V. Ruuskanen, and 
K. Tuominen, Nucl. Phys. B570, 379 (2000). 

[17] R. Paatelainen, H. Holopainen, and K. J. Eskola [in 
progress] . 

[18] J. P. Boris and D. L. Book, J. Comput. Phys. All, 38 
(1973). 

[19] S. T. Zalesak, J. Comput. Phys. A31, 335 (1979). 
[20] M. Laine and Y. Schroder, Phys. Rev. D73, 085009 
(2006). 

[21] S. S. Adler et al, Phys. Rev. C69, 034909 (2004). 
[22] K. Aamodt et al. [ALICE Collaboration], Phys. Rev. 
Lett. 105, 252301 (2010). 



[23] J. I. Kapusta, P. Lichard and D. Seibert, Phys. Rev. D 
44, 2774 (1991) [Erratum-ibid. D 47, 4171 (1993)] [Phys. 
Rev. D 47, 4171 (1993)]. 

[24] H. Holopainen, S. S. Rasanen, and K. J. Eskola, Phys. 
Rev. C 84, 064903 (2011); M. Dion, C. Gale, S. Jeon, 
J. -F. Paquet, B. Schenke and C. Young, J. Phys. G 38, 
124138 (2011). 

[25] R. Chatterjee and D. K. Srivastava, Phys. Rev. C 79, 
021901(R) (2009); R. Chatterjee and D. K. Srivastava, 
Nucl. Phys. A830, 503c (2009). 

[26] K. J. Eskola, K. Kajantie and K. Tuominen, Phys. Lett. 
B 497, 39 (2001). 

[27] P. F. Kolb, U. W. Heinz, P. Huovinen, K. J. Eskola and 
K. Tuominen, Nucl. Phys. A 696, 197 (2001). 

[28] K. J. Eskola, P. V. Ruuskanen, S. S. Rasanen and 
K. Tuominen, Nucl. Phys. A 696, 715 (2001). 

[29] T. Renk, H. Holopainen, R. Paatelainen, and K. J. Es- 
kola; Phys. Rev. C 84, 014906 (2011). 



